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ABSTRACT 



Aq experimeD,tal aad theoretical analysis was 
performed on low e-nsrgy impact loading of thin (.040 
inches) graphite-e<pox y plates. Six-inch sguara plates 
were subjected to dynanic impact^ (below the ballistic 
range) and static loading. The plate static strain 
energy and dynamic impact energy for failure were 
equal and constant. Theoretical analysis was 
performed using both exact and finite element methods. 
Small deflection theory was assumed and found to be 
inapplicable; the plates behaved in a manner 
indicating that a large deflection theoretical model 
would be more appropriate. 
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I. 



INTRODUCTION 



Composite materials have a tremendoas potential in the 
structural field. Although their future value in the 
aircraft industry is exceptional, there 5 till remain many 
unknouns. The response of graphite-epoxy composite plates 
to low energy impact loading is one of these unknowns. 

Low energy impacts are a common occurrence on 
operational aircraft. A mechanic dropping a wrench, or a 
stone thrown by a tire are typical examples. It is 
important to know how composite panels react to these 
impacts in order to be able to confidently predict a damage 
threshold. An accurate model for damage would be helpful in 
both the operational and design fields. Operationally, the 
ability to assess damage to a panel which is known to have 
been hit, but which has no visual manifestations of damage, 
will produce better maintained and safer airoraft. The 
designer will also be better able to produce aa aircraft 
capable of operation throughout its design environment. 

A series of tests, both static and dynamic, were 
performed to analyze the response of a simply supported 
plate to impact loadj.ng. The plate was a symmetrically 
laid-up graphite-epoxy plate six inches sguare. The dynamic 
tests used a variable-<-mas3 steel penetrator of spherical 
shape, at impact speeds of ten, fifteen, and twenty feet pet 
second. Each plate was impacted by incrementally increasing 
masses until failure occurred. 

There were three basic types of failure expected. These 
included: static shear-out, where the penetrator sheared 
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through the panel causing only local damage; flexural, where 
the failure was caused primarily by bending aad membrane 
stresses; and ballistic, where complete penetration 
occurred. It was assunei that wave stresses would have 
little influence, and that flexural failure would 
predominate. Saall deflectioas were iaitiaily assumed and 
membrane stresses were neglected. 

The Kirchoff Hypothesis for plates was adopted: 
straight lines normal to an undeformed middle surface were 
assumed to remain straight and normal to the middle surface 
in the deformed state. Further, it was assumed that normals 
to the middle surface do not change length. 

The investigation proceeded in the following manner: 
first, the constitutive properties of the plate material 
were determined both analytically and experimentally, then 
the experimental apparatus was sat up and calibrated to give- 
the desired range of impact velocities. k number of 
non-destructive technigues for inspection of damaged plates 
were considered and tested. Finally, a series o£ static and 
impact tests was run, and the resulting damage was evaluated 
using both nondestructive and destructive test methods, 
including scanning-elect con- microscopy. 

Concurrently a theoretical investigation was undertakec 
to see whether simple linear theory could be of use in the 
analysis of the plate impact problem. Static aad dynamic 
solutions were sought using SAP IV, a finite-element 
structural analysis program. Initially, known closed form 
solutions were compared with SAP IV solutions for static and 
dynamic loadings. Then specimen constitutive properties and 
experimental loading parameters were entered in tie program, 
and solutions were compared with experimental data. 

The procedures used and the results gained ace detailed 



13 



in 



the following 



sections. 
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II 



ZDfJSTirUnVE PROPERTIES 



A. MATERIAL PROPERTIES 



The plates used were roDstructel frem 12 iash wide rolls 
of UCC Thornel 330 (uatwisted) Sraptite Fibers, in a 
Rigidite 5208 matrix, supplied as a B-stage prepreg by the 
NARMCO Materials Division of the Celanese Corporation. 

Twelve-inch sguare panels ware constructed following the 
manufacturers recommended procedures, and those of Linnandec 
Tfief. 1]. These proceiures were designed to give a 55i5 
fiber content by volume. The lay-up was eight lanina in a 
symmetrical configuration as indicated by Fig. 1. The 
standard laminate orientation code for this lay-up [Ref, 2], 
was (0,±45,0) . After the proper cure, the panels were 



s 

guartered to give 



six-inch 



square plates. 



Some of the lamina properties predicted by NARMCD 
Materials included, in pounds per square inch: 

Flexural modulus (longitudinal) =22 , 3 B6 
Flexural modulus (transverse) = 1.6 E6 

Tensile modulus of coupon sample (logitudinal) =23. 5 E6 
Tensile modulus for coupon sample (transverse) =1 . 47 E6 
Flexural ultimate stress ( longitudinal) =295 E3 
Flexural ultimate stress (transverse) = 11.5 E3 
Tensile ultimate stress for coupon sample 
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Figure 1 - 0SIENTRrr3N OF LRMINR ON R LAMINRfSD PLATE* 
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(longitudinal) =2 1 5 E3 

rensile ultimate stress £dc coupon sample 
(transverse) =6 . 75 E3 

Compression ultimate stress (longitudinal) = 193 S3 
Interlaminar shear ultimate= 17.8 S3 

The published properties did not iaclude either of the 
Poisson ratios or the shear modulus. The nominal thickness 
of a laid-up plate was .3*3 inches. 

B. THEORETICAL C0NSTITUrE7E PROPERTIES 

Jones [Ref. 3] developed equations for the constitutive 
properties of a laminate from known lamina properties and 
laminate lay-up. 

Beginning with the known properties the constitutive 
properties of each lamina in relation to the Laminate was 
determined. Lamina with Eiber direction coincident with the 
axis of the plate; and plane stress: 




For laminas that were not on the principle plate axis the 2 
matrix was modified with a suitable coordinate 
transformation. This applied to the t 45 degree Laminas: 

= Q» COS‘*e -t- a(Q,i + 2 Ofc4.\Sl M^C0S*9 + OzzS \ 

= fQa+Q22-HQfofe)S\N"ecofe+Q,2(SlW‘'0^CO5‘^©) 

Qa = QuSUvi^'G + 2(Qii-^ 2Q4.fe) 5IM*0CD5^G 4-OiiCOs'*e 

Qife = (Qu-Qii-2Qa.^51MecoS^0 + ('O,rO22-^2Qu')5»t^ecose 

SL ^ (On -2Qfc<.^5lM?0CCS0+ (Qit-0u+ZQfefe^51M0CO5’© 
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Qfeu = CQ,,+Qa -ZQu-2Qu'^SlN^eCDSb + a4f5llvf©-^C0S'^©') 

Foe small strains and linear elasticity, utilizing 

classical laminate theory, the stresses in any Layer were 

expressed in terms of middle surface strains and curvatures. 

By integrating the stress in each lamina over the thickness 
of the laminate: 



Ny ) ♦ 07 dz 
X X 



and 



(Xzdz 
^2 

the following relations were ootained; 



M 



X 



^n] = [a J ^ [b J (x] 

- [b ] {e-] + [d ] ^x} 



where for extensional stiffness; 
for coupling stiffness: 



B. 



{Q„ ).. (z!. -Z^ 



Ij ii K K 

for bending stiffness 

D;; = W > 



K-l 



IJ 



K* 






The plates as constructed were symmetrical, 
the B, or coupling, matrix to zero. 



This forced 



To facilitate calculation it was desirable to express 
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the constitutive properties in the form of engineering 
constants of an octhotropic material. The finite element 
program used woild accept data in this form only. 

The Design Manual [Raf. 2] supplies a means to 
approximate the ortiotcopic constitutive properties 

utilizing the k, or extensional stiffness, matrix; 

2 

Exx = (All A22-A12)/ (THICKNESS)A 22 
^yy “ ^xx ^22^ A|| 

Gxy -A33/ THICKNESS 
yx “ A |2 /A|| 

^xy ’ A22/A12 

A simple computer program was designed to ta^ce both the 
lamina constitutive properties and laminate construction as 
input. The program, in Basic Language, outputs the 
following: the individual lamina stress-strain properties 

in relation to the laminate, the extensional matrix, the 
bending matrix, the orthotropic constitutive properties, and 
the orthotropic stress-strain relationship. This program is 
Appendix A. 



C. EXPEfildENTAL CONSTITUTIVE PROPERTIES 



To verify the calculated constitutive properties, a 
series of tensile tests was performed. 

Three specimens were prepared in accordance with the 
design manual [Ref. 2]. Uaey ware nine inches long and one 
inch wide. The first specimen was cut from a twelve-inch 
panel with the zero degree fibers running parallel to the 
long dimension. The second was cut so the caro degree 
fibers ran at a right angle to the long dimension. The last 
specimen was cut so one of the forty-five degree lamina's 



19 



fibers were parallel to taa long dimension, and the zero 
degree fibers were fort/-five degrees off the Long specimen 
axis. One-and-a-half-inci long and one-inch wide shoulders 
were glued to both sides of the outer one-and-a-ialf inches 
of the specimen in order to distribute the local stresses 
caused by mounting in the test machine. The edges, where 
the shoulders met the specimen, were beveled to less than 30 
degrees. 

3BH Electronics SR-'4 strain gage rosettes were installed 
on each specimen near the center. These rosettes were type 
FAER-25R- 1 2S6 with 120 oin resistance and 2.04 gage factor. 
They were wired into a standard Wheatstone bridge strain 
measuring circuit. 

The three specimens were mounted in a RIE3LE testing- 
machine. As they wars slowly loaded to failure, strain 
measurements ware incrementally taken. Figures 2, 3, and 4 
are plots of the test data. 

The orthotropic constitutive properties can be 
determined from the graphs. From this data, the Design 
Manual [Ref. 2] gives a simple calculation procedure to 
determine the shear modulus: 



U*(EL-^ Et-I-2S)ltEy)/8(I-Vlt''tl) 




2 ^ Ex 



From the graphs: 



El* 11.65 E 6 
Ei*» 3.1 E6 
Ey - 61.6 E 6 
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Figure 2 - TENSILE TEST DF ZEKD nE3R-BE-cS.PESIMEN 
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Figure 3 - TENSILE TEST DF 90 DESREE SPESIHEN 
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To verify Poissons ratio: 

IlT-'^JL - « .14 

The shear modulus calculation gave: 

Glt= 2.0 5 E 6 

Summary of constitutive properties from tests: 



El= 11.6 E6 




E-r= 3.1 E6 


\t'-52 

V \ 


Glt= 2.1 E6 


Mtu-.I4 



Using this data the stress-strain relationships will be: 




D. COMPARISON OF CONSTITUTIVE PROPERTIES 

Since Poisson's ratios and tae shear modulus for an 
individual lamina were lot given, these properties were 
estimated and entered into the constitutive properties 
computer program. Using the orthotropic properties obtained 
from the tests, the program was used to iterate for the 
unknowm lamina properties. The results were disappointing. 
The calculated modulus in the 5f direction was eight percent 
high, the shear modulus was forty percent high, and the 
Poisson's ratios were thirty five percent high when the X 
direction calculated modulus was forced to equal the test 
data. 



A survey of the Design Manual ‘8ef.2] indicated that the 
constitutive properties obtained experimentally were 
comparable to common, medium strength graphite epoxy places 
of similar construction. The orthotropic values obtained 
from the test data were used unless otherwise noted. 
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E. 



FIBER-EPOXY VOLUME FHAOriON 



Two independent tests were conducted to verify the fiber 
and epoxy volume fractions. One method was x chemical 
technique, the other was a quantitative microscooic method. 

Handley and Cross TBef.i] outlined the chemical method 
used. Two fritted crucibles and two samples were dried in 
an oven at 180 degrees centigrade for thirty minutes. The 
crucible and coupon were then weighed. Tne coupon was 
placed in 80-90 degree centigrade, concentrated nitric acid 
for fifteen to twenty minutes. The hot nitric acid 
dissolved or "digested" tne epoxy matrix. After the epoxy 
was dissolved , the reniining fibec and acid sol ut ion was 



rinsed with a 


vacuum assist in 


wate 


r and 


acetone 


The 


crucible with 


the fibers was dried 


in 


an oven 


at 120 


degrees 


centigrade for one to two hours. 


and 


then 


weighed 


The 


results were: 












Sample One 








Sample 


Two 


1.0612 


specimen weight 


in 


grams 




1. 3525 


14.1733 


specimen and crucible 


(before) 




1 4. 2385 


13.9225 


specimen and cruciole 


(after) 


1 3. 9307 


.2508 


epoxy weight in 


gcaais 




. 3378 


23.6% 


epoxy by 


weight 




24.8% 


1.265 


specific density 


of 


epoxy 




1.255 


1.9-2. 3 


specific density 


of 


graphit 


g. 


1 .9-2. 3 


54-68% 


fiber bv 




volume 




63-67% 



This compares favorably with the desired 65%. 

A quantitative microscopic analysis was done to further 
verify the fiber/epoxy mixture. 

A section was taken from a plate using a diamond cut-off 
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vt(k 



saw, and cold mounted in plastic. tt was than polisnei 
using grinding baits of 3D and 320 grit, emory paper of 0 
and 3/0 grit, and polishing cloths with aluminum-oxide 
slurries of 15 and .05 micrometer grit. The microstr ucture 
showed sufficient resolution with just polishing, so etching 
was not used. A nitric acid etch would have selectively 
attacked the epoxy, it desired. The cold mounting process 
caused slight sample degradation. The plastic attacked the 
epoxy in the outermost layer. The attack was very slight 
and confined to dissolviag the extreme surface layer of 
epoxy. If there was a well developed crack, the plastic 
would also fill these with a corresponding light attack on 
the surrounding matrix material. The plastic mounting 
material was also significantly softer than the graphite 
epoxy plate making it difficult to obtain sguare edges on 
the sample. 

Figure 5 details the results of the guantitative 
analysis. The illustration is the 3DDX nicrophotograph used 
for the analysis, which is a picture of the middle two 
layers of a plate. The data indicates a fractional length 
of Line within the finer as .557. This eguates to the 
volume fraction giving a 35.7^ fiber volume, which concurs 
with previous analysis. 

The average dimension of the fibers was also determined 
from the microphotogcaph using standard guantitative 
microscopy tecnnigues. The average dimension was 5.4 
microns. Since the specimen was cut at a right angle to the 
principle fiber orientation, this dimension gives a close 
approximation of the average fiber diameter. 
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Figure 5 - MICROSCDPi: ANALYSIS 3F FIBER-VDLU'IE CONIENr 




III. 



EXPEftlMENIM. APPRAiaS AMD P2SLIMINARI TESTS 



A. EXPERIMENTAL SET-UP 



1 • PeEetratgr 

Two variable mass psaetratocs wacs constcucted and 
are illustrated in Fig. 5. The primary penetrator was the 
lighter, slimmer one. It was constructed with a plastic 
pipe body sized to fit in the launcher tube. The mass could 
be adjusted from approxinently .314 to .393 slugs. This was 
accomplished by varying the weight from seven ounces to 
three pounds with lead shot. The leading point of the 
penetrator is half-spherical shaped and made of steel. The 
larger all steel penetrator had a mass which could be varied 
starting at .085 slugs. This penetrator had the capacity to 
change both the leading point's shape and its construction 
material. The second penetrator was not needed in the 
dynamic tests. 



2. Plat e Fix-tur e 

A massive steel fixture was built for both the 
static and dynamic tests. This fixture is illustrated in 
Fig. 7. It was carefully machined to have a flat test 
surface and to hav*e the top and bottom surfaces parallel. 
The fixture was made from heavy steel channels with a 
one-naif -inch thick steel base. 
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Figars 5 - PENErRAIDRS 
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Figure 7 - PLAIE lESriNG FIXTURE 



30 




i 



Tha square test supports measured four- a ad-one- half 
inches on the inside. Using a six-inch plata, aach plate 
edge over-lapped the support by three quarters of an inch. 
The fixture also had an additional frame and a series of 
clamps to secure all edges of the plates to the fixtures; 
this frame would prevent all edge rotation. Plates could be 
mounted and tested in either the simply supported or clamped 
edge modes. Large holes were drilled in the sides of the 
fixture to prevent any pressure build-ups auriag dynamic 
displacements. 



3 . Launcher 



The desired impact velocity was achieved by simply 
dropping the penetrator from an appropriate height. To 
insure that the penetrator would hit the desired spot on the 
plate, and to provide for repeatable, easily adjustable drop 
heights, the launcher snown in Fig. 3 was constructed. The 
launcher could be varied in height above the plate via a 
screw crank for fins velocity adjustments. Both penetrators 
were just slightly smaller in diameter than the inside 
dimension of the launcher tube, which was used to stabilize 
the penetrator's flight. Launch was initiated by a switch 
controlled solenoid waich removed a bar holding the 
penetrator in place. 

4 . Frame 

A frame was constructed out of steel scaffolding, 
and is illustrated in Fig. 9. The frame vas plumbed 
vertical and constructed on a heavy-load rated concrete 
floor. The launcher assembly was moved up and down the 
frame by a pulley and rope arrangement to obtain course 
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velocity 



ad j ustmea ts. 
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Figure 3 - PENSIiJArOR LAJNCHER 
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Figure 9 



LAUNCHSR FRAME 
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5 . 



Static Tests 

Static tests wsce coaductei on a P.IEiLE testing 
machine using the base tixtura, plate and appropriate 
penetrator. Deflections were measured with a dial 
displacement indicator from the cross am of the testing 
machine to the working surface. 



B. CALIBRATION 



The launcher-frame coabinatLon was calibrated for both 
impactors at three different speeds. Ten feet per second 
was the lowest practical speed obtainable due to launch tube 
interference. Twenty feet per second was the highest speed 
desired in order to keep in the flexural failure mode. 
Fifteen feet per second was also used. 

The calibration was done electronically . Two light 
sensitive photo-transistors and their power supplies were 
placed two inches apart, bracketing tie plate fixture's 
position relative to the launcher. Two ligat beams ware 
focused across the penetrator flight path onto the 
transistors, and an electronic timer was placed oetween the 
transistors. The wooden frame supporting this apparatus was 
painted flat black to prevent spurious imputs. Speeds of 
ten, fifteen, and twenty feet per second were calibrated 
within two percent. 
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c. 



IMPACT DURATIjN 



The time and strain-amplitude tiistory Df a typisal 
impact was obtained to assist in the theoretical 
calculations . A graphite epoxy plate similar to the tested 
plates was fitted with a strain gage near the center of a 
face. The strain gage was connected to a Wheatstone bridge 
and power supply. An oscilloscope was connected across the 
bridge and set to trigger on impact. The scope was set for 
single sweep mode, and had a camera mounted to record the 
trace. 

Figure 10 is a scope picture taken using the larger 
impactor at its lightest weight. The spaed of impact was 
approximately ten feet per second. The oscilloscope time 
scale for this photograph was .5 miliseconds per centimeter. 
The amplitude was not calibrated and gives relative 
information only. The strain in the plate reached its peak 
value at one millisecond according to this photograph. 
After the impact, examination of the plate determined that 
the strain gage leads had separated from the gage leaving 
the validity of the results in doubt. The test did give 
several useful pieces of information. It gave an indication 
of the impact’s initial siape and rise time, and, since wave 
stresses would have reacied maximum values in much quicker 
time ranges (microseconds) , the test confirmed the 
assumption that a flexural or membrane mechanism would be 
the primary loading mode. 
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Figura ID - TYPICAL IMPACT TIME HI3T3SY 
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D. 



NON-DESTRUCIIVE TESTIMS 



It was necessary to examine each plate after every 
impact loading to determine if any damage had occarred. The 
tests had to be quick, inexpensive and non-destructive since 
many tests would be required during the course of the 
experiments. Several methods of non-destructive testing 
were considered, including: ultrasonic sound. X-rays, and a 
new liquid crystal process. 

1 • Ultrasonic Testing 

An Automacion Industries Inc., Sperry Division, 
Model QM-771 Ref lectoscope Ultrasonic Testing )lachina was 
examined for its ability to meet tne requirements. The 
machine was equipped with 2.25 MHz and 5.0 MHz zero degree 
transducers. It also had 2.25 MHz transducers with sound 
input angles of 45, 50, and 90 degrees. All the 
transducers, including various combinations, were used on 
two different graphite epoxy panels with known flaws. One 
panel had the standard surface finish, and tne other was 
specifically manufactured to give a very smooth surface 
finish. Both the shear and through modes were attempted, 
along with the different angle input modes. 

The panels with normal surface roughness, from the 
bleeder plies used in manufacture, were difficult to test. 
The sound return from both the front (input), and back 
surfaces was lost in the clutter caused Oy the roughness. 

The plate with a smooth surface allowed a reasonable 
input and back surface echo, but the plate was so thin tnat 
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these ceflectioas occupied the entire useful range of the 
echo returns. Known flaws could be detected in the clutter 
using the highest fraguenc/, but it was difficult to get any 
continuous, consistent returns. This was true even with 
known, large flaws. 

Other modes of the machine were tested with negative 

results . 



For a plate as thin as tais, a higher frequency 
transducer might have given useful results, but for the 
purposes of this investigation ultrasonic tasting was 
discarded. 

2. X-Ray Testing 

A portable X-Ray dachine was triad as a means of 
non-destructive testing. A General Electric LS-40 Portable 
Industrial X-Ray Jnit was set-up to examine a graphite epoxy 
panel with known flaws. Polaroid film was exposed at 
various kilovolt and time settings to determine the 
procedure's usefulness. Even at the lowest setting of 70 
kilovolts and 10 seconds, the plate did not absorb enough 
energy to produce an inage. The pictures obtained had no 
contrast and no detail. k machine capable of even smaller 
energy output would probaoly still have difficulty obtaining 
any contrast due to the light atomic weight of the carbon. 
This means of non-destructive test was abandoned. 

3* Cr ystal T§.3ting 

Shaum [Ref. 5] developed a method for using 
temperature sensitive liqaid crystals for damage evaluation 
in graphite epoxy laminates. This method was used with 
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exceileat results. 



After prepariag plates for testing, they were 
painted flat black on aie side. Liguid crystals, with 
temperature sea si tivit y oi 30±1 degrees centigrade, were 
mixed with water and a snail guancity of wetting compound 
(liquid detergent) . This mixture was then spray painted on 
the plate to be tested. 

After each impact during the dynamic testing, the 
plate was examined for damage. This was accomplished by 
selectively heating all, or just a portion, of the panel, 
and watching the heat transfer cause the liquid crystals to 
change color. Any damaged area would experience a different 
heat transfer rate and thus would change color at a 
different time from an undamaged part. This method was 
sensitive enough to show even slight compressive 
indentations in the surface prior to their becoming visually 
or tactually evident. 
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17. FE.UIE ELSHSNr ANALYSIS 



A. FINITE ELEMENI MODEL 



SAP IV, a structucal analysis program for static and 
dynamic response of linear systems by Bathe, /filson and 
Peterson [Ref. 5] was used as an analytical tool. 

Two finite element nadels of a plate were made for 
analysis. The first model was a coarse model coasisting of 
only 16 square plate elements, 1.125 inches on a side, for a 
4.5 inch square plate. This model was used for initial 
"coarse" investigation at a small cost in computer time. 
The second model was a fine model with 324 plate elements, 
.25 inches square, and used considerable computer time. 

The results produced by SAP IV using these models were 
compared with known, exact solutions for aluminum plates. 
Timoshenko and Woinowsky- Krieger [Ref. 7] solved the simply 
supported, concentrated load problem, which was used as the 
reference solution. Figure 11 shows a plot of the 
deflection versus position on a plate (from an edge to the 
center keeping I constaati under a 103 pound load. The fine 
finite element model is . -t % stiffer than the exact solution, 
where the coarse model is 5.5% stiffer than the exact 
solution . 

The dynamic analysis section of SAP was addressed to 
compare its modal freguency with an exact solution. 
Meirovitch [Ref. 3] gives an exact solution for the 
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vibration frequencies of a recfcangulac plate: 



FREQUENCY** 







m,n ** 1,2,... 



Ihe coarse finite element model results were within .5% of 
the exact solution. Since this agreement was considered 
more than satisfactory, cae fine model was not rua. 



% 
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EXfla RND FINITE aEHENT (FINE) HnLUTlDNH 
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igure 11 - EXACT SCLOTION COMP&RED WITH FINIfE ELErtENT 
SOLUTIONS FOR A CONCENTRATED LOAD ON AN ALUHiyiUH PLATE 



I 




B. STRESS DISTRIBUTION 



In order to limit tte required amount of ooiputer data 
reduction by determining where initial failure could be 
expected, an analysis was performed using 5AP IV to 
determine the stress distribution on a typioal plate. 
Figure 12 illustrates the results of this study on a simply 



supported plate 


with 


a concentrated 


load. r 


he 


stress 


concentration was 


the 


highest at the 


center of 


the 


plate. 


and this indicated 


that 


initial failure 


analysis 


CO 


uld be 



Limited 

plate. 



to the immediate vicinity of the center of the 
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DISTANCE FROM EDGE TO CENTER OF PLATE (inches) 
Figure 12 - DI!i3aAM 3F SrSESS DISTRIBOriON 
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FAILURE CRIIERIDN 



rests made b/ McQuilLsQ/ Lloceas, aai Gause '3ef. 9] on 
composite beams demonstcated a flexacal failure mode. Thus, 
for this investigation the assumption was made that the 
flexural failure mode woiLd predominate. Since sach lamina 
was significantly weaker in the transsrarse direction than in 
the direction paralLal to the fibers, and since the 
outermost layer would have the highest stress leval, it was 
expected that the preliminary failure would be in the 
outermost layer ia the transverse direction. After this 
failure, the load normally carried would be transferred to 
the next layer. If this lamina could not carry tne load in 
its transi^erse direction, failure of the plate was assumed. 
Bach lamina’s longitudinal strength was so large that no 
failure was allowed in that direction. 

Returning to Jones ' Xef . 3], the strains are: 



Since it is assumed that tie length of the middle surface of 
the plate is constant, the first partial terns can oe 
dropped . 

Using, 



the middle surface curvatures are found. Having these, and 
the stiffness matrix for each lamina, the lamina stresses 






are fou’^'^ fi-nm* 
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Knowing the stresses ani the ultimate failure stress, 
failure can be predicted. 

A computer program was written in Basis Language to 
evaluate failure. The input consistei of the moments per 
unit length in the X, X, and XX iirections. Ehe program 
started by computing the stresses in the outer two layers, 
as outlined above, using appropriate coordinate 
transformations. It then set the transverse and shear 
modulus for the outer lamina to zero and computed the 
stresses on the next layer. Failure was assume! to occur 
when the transverse ultimate failure stress was exceeded in 
both the outer and next layer. The program is included as 
Appendix 3. 
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V. 



STATIC TEST RESULTS 



Seven static tests wece performed, six involving loading 
to failure. Three tests used the larger penetrator and 
simply supported boundaries. Two used the smaller 
penetrator and were also simply sapported. Dae test was 
conducted using fixed boundary conditions and the small 
penetrator. The last plate was not loaded until failure; 
instead, it was partially loaded and then used for 
microscopic examination. 

Figures 13 and 14 show the results of all static tests. 
The static failure energy, as outlined by ScQuilien, 
Llorens, and Cause 'Ref. 9] for beams, was one half the load 
times the displacement. 

All failures, with one exception, involved the same 
failure mode, and Figure 15 and 15 are typical examples of 
this mode. The only exception was the test #ith fixed 
edges, and Fig. 17 illustrates its mode. 

During the course of the static tests it was apparent 
that a "yield" point was experienced. This point occurred 
at approximately half the ultimate load when the panel would 
experience a decrease in load with no increase in 
displacement at the center of the plate. figure 14 
indicates that this was common to both penetrators and 
boundary conditions. Conseguently , a plate was loaded to a 
point just beyond this "yield" point, and then sectioned and 
prepared for microscopic examination. As the plate was 
loaded a large deflection occurred in the plate which was 
very similar in size and snape to the standard failure. 
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Figure 18 indicates this ieflection area, and shDrfs were the 
sections (marked 5) were taken. 
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Figurs 13 - TABLE OF STATIC TEST RESULTS 
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Figura 14 - 3 RAPH OF STiirrc TEST RESULTS 
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Figure 15 - SKETCH DF A TYPECAL STAriC OR DYNAMIC FAILURE 




Figure 16 - E333 VIEW DP k liPICAL FAILaRE 
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Figure 17 - STATIZ 3HBAR-0Ur 
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VI. 



DlfSAMIC TEST RESULTS 



Figure 19 simmarizes the results of the iyniinic tests. 
Each plate was iapactad repeatedly by tae panetrator, at 
constant velocity but with progressively increasing mess 
until failure occurred. Failure was evidenced oy a loud 
cracking noise, virtually no penetrator rebound, and gross 
delamination. The energy reguicad for failure is indicated, 
and is one-half the mass times the velocity sguaced. After 
each impact the plates were examined for damage using the 
liquid crystal inspection technique. In all of these tests 
the smaller size impactor was used. 



Each 


failure. 


with one exception, shows 


d 


the 


same 


failure 


mode as the 


static case shown in Fig. 


13 


and 


16. 


The one 


exception 


occurred at 20 feet per 


s e 


con d. 


and 


approached ballist 


ic in nature. Figure 20 


shows 


this 


failure 


with its 


almosc complete penetration. 




splaying 



of the delaminated layers on the back (non-impacti side. 

Early in each test, damage was observed in the plates. 
This damage was localized in the immediate area of impact. 
When it started it was oarely discecnible using the liquid 
crystal technique, but as the tests progressed it became 
easier to see. After many impacts, the damaged area was 
delaminated, and this could be determined by careful tactile 
examinations. It did not appear to have any effect on the 
the ultimate failure energy. This '•bruised" area was 
subjected to a microscopic examination. It was section 1A2 
in Fig. 18 and is shown enlarged in subsequent figures. 
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TEST SPEED NUMBER OF PENETRATOR ULTIMATE 

(ft/sec) IMPACTS UNTIL WEIGHT IMPACT ENERGY 

FAILURE (lbs) (ft-lbs) 
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Figure 19 - TftBLE OF DiTNA.IIZ TEST RESULTS 



VII. 



jrSCOSSIDN DF RE5ULTS 



There were three striking phsnoaena encountered during 
the coarse of the tests. The first was that the strain 
energy for failure in the static tests and the impact energy 
for failure in the dynamic tests were egual and constant. 
The second was the initiiL dynamic damage that occurred well 
before the ultimate failure. The third was the large 
deflections encountered. 

as the tests and analyses progressed, the large 
deflections had an important influence on the correlation of 
experimental and theoretical results. The plates were .J40 
inches thick and static failure occurred at plate center 
deflections of approximately .35 inches. This means 
deflections were eight to nine times larger than the 
thickness of the plates. 

An analysis was attempted to correlate the actual 
results with theoretical solutions. Five exact and finite 
element solutions were ootained for the static case and an 
intensive structural dynamic response analysis was attempted 
for the dynamic cases. These investigations were thwarted 
by inadequacies in the basic assumption of small 
deflections . 

The plates all experienced a common failure mode for 
both the static and dynamic tests. microscopy was utilized 
to examine this ultimate failure surface and mechanism. 
Further areas addressed with microscopy included: a 
pre-failure dynamic damage zone, or "bruise," tnat developed 
on each panel; a "yield* point that occurred during static 
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loading at approximately oaa half the ultimate load; and the 
fiber failure mechanism. 



A. STRAIN RATE DEPENDENZE 



As indicated in Fig. 14, in ona large-panatr a tor tast 
the specimen seamad to ha significantly stiffac than the 
other plates, all of which showed guita similar 
load-deflection curves, cagardless of paaatrator size. The 
only parameter which differed in this "stiffer" tast was the 
rate of loading, which was higher than that jsed in the 
other tests. No further gualitativa analysis was performed, 
but the rate of loading was kept as low as possible. 
Further tests then yialdad consistent, repeatable results. 



B. FAILJRE MODES 

All failures, both static and dynamic, showed the same 
failure modes [Fig. 15, and 16]. There were only two 
exceptions. The first was during testing of the fixed plate 
with the smaller oenetrator. This failure was illustrated 
in Fig. 17, and was a form of saear. Although the plate was 
stiffer than the simply supported plates, it responded in a 
similar manner until the "yield" point was reacaed. After 
that point, the penetrator virtually sheared tarough the 
plate. The second exception was one test at 20 feet per 

second. The failure mode approached ballistic in nature. 
Figure 20 illustrates the almost complete penetration, and 
the splayed lamina around the exit hole. 

The choice of 20 feet per second as a limit for the 
dynamic tests was justified. Since one of the four tests 
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conducted expecisnced this failure mole, this speed range 
was in the transition fcoa flexural failure node to a 
ballistic failure mode. 

rhe size of the plate was sufficient to keep the place 
in the flexural failure, instead of shear failure, mode 
during static tests. All tests in the simply supported mode 
were flexural in- nature. 



iliS£ 2 scopic Exaniaation of Oltimate Failure Area 



Section 1A3 [Fig. 18] was cat and mounted for light 
microscope examination. Figure 21 shows the region outside 
of the gross delamination area typical of failure. This 
section appears relatively intact with vary minor 
delamina tion . There ace several small lamina fractures in 
the outer two layers running with the fiber directions in 
the matrix, and these minor cracks are associated with a 
slight delamination. 

Figures 22, 23, and 24 show the lamina in the gross 
delamination regio n. These show extensive de lamina tion, and 
numerous large cracks in tne matrix, between fibers, 
penetrating entire laminas. Each delaminatioa between 
laminas appears to have at least one crack penetrating the 
lamina directly above or below. Figure 24 illustrates this, 
with delaminations and associated fractures of tie adjacent 
laminas . 



Totally flexural failure damage should have been 
concentrated in the outermost layers where tne largest 
stresses occurred. Figures 22, 23, and 24 all demonstrate a 
significant failure of the middle two layers in their 
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Figure 21 
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Figure 22 



LISHI i1ICR3320PE, FAILJRE SURFACE 
SECTION 1A3, 130X 



EXAMINATION, 
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Figure 23 



LI3HI HICaDSZDPE, FAILtJSE SQRFA3E 



SEZnON 1A3, 53X 



SKAMINArrON, 
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Figure 24 



LI5HT MI3R33:0PE, FAILaRS SURFACE 
SECTION 1A3, 2D0X 



EXAMINATION, 
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transverse direction. rtie assumption that the middle 
surface of the plate iid not experience extension, 
therefore, is invalid. Tais would indicate that cheoretical 
analysis based on small deflection theory is inappropriate. 

In Figures 21, 22, and 24 the damage is most intense 
on only one side of the plate. The outer two layers on the 
side opposite the applied load appear almost intact. A 
quantitative examination of all failed plates showed an even 
distribution of damaged outer layers. Some plates showed 
gross lamina fracture in ooth outer layers. Dtier panels 
had damage on only the loaded side and vice-versa. Thus, 
the two intact layers are not considered significant. 

A section was cut for the scanning electron 
microscope (SEM) evaluation as indicated in Fig. 13. Since 
this was an area of gross delamination, the top (or loaded 
side) two layers separated during sectioning. These layers 
were prepared for evaluation. 

Figure 25 shows a typical fracture, or crack, in the 
epoxy matrix running tnrough the outer layer. Figure 25 
demonstrates a similar fracture, further illustrating the 
complete penetration of the lamina. The entire outer layer, 
a zero degree layer, was extensively cracked. Figures 27, 
28, 29, and 30 show typical fractures in the first 45 degree 
layer . 



Figures 27, 29, and 30 have remnants of the third 
inner layer (45 degrees) which are evident as epoxy and some 
fiber fragments running at right angles to the original 
layer. This was typical of all the gross de lamination 
surfaces . 
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Figure 25 



5Ed, rypEZAL FRACrJRE IN OUTER LATER, 523X 
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Figure 27 - SEM, TYPE:&L FRACTURE IN 5EC0ND LAYER, 1200X 
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Figure 28 



SEM, TYPEZM FSACrORS IH 3ECOSD LAYER, 2200X 
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Figure 29 



SEM, TYPIZM FRACTURE IN SECOND LAYER, 



1200X 
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Of particular iatscest is the contiuuity of the 
fibers in the layers. ILL inter-L aninar failures occurred 
at the epoxy-fiber joiats with no discerniole fiber 
breakage. A quantitative analysis, presented in Fig. 31, 
indicates that slightly over ten percent of the inner lamina 
epoxy-fiber joints have broken. These failures are long, 
parallel to the fibers, continuous, and appear to penetrate 
completely. This would indicate an inability to carry any 
load in the transverse direction while still naintaining 
strength in the fiber direction. This would be consistent 
with the low transverse, but vary high Longitudinal, 
strength and failure char acter ist ics. The transfer of this 

the next lamina would cause large 
forces, and this could account for the 
The observation that all deL amin ations 
associated matrix fracture, or crack. 



transverse load to 
inter-laminar shear 
gross delaminations, 
appear to have an 



further supports this theory. 
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DF FAILURE 3F FIBER-EPDXI JOINTS 




I 



2* lll££2§£2Ek£ ExasLaacion of Static HX,i§liH PoiSi 



During static loading of the plates, there ^as a 
pronounced "yield" point at appcoxinately half the ultimate 
load. At this point the panel would experience a decreased 
load with no increase in displacement, accompanied by loud 
cracking noises. Figure 14 gives plots of these points. A 
plate was loaded just past this "yield" point, and then 
sectioned (5B) as snowa in Fig. 13. Using the light 
microscope. Fig. 32 shows an area where all eight lamina ace 
delaminated. In Fig. 33 there are just a few delaminations, 
but both pictures show extensive damage to the underside 
(side opposite the load} zero degree layer. 

Figures 34, 35, aad 36 all show the same damage to 
the bottom layer along with assorted other delaminations and 
cracks. Figures 34 and 35 have matrix fractures associated 
with delaminations. In the microphotographs, the damaged 
outer layer was probably accentuated by the reaction of the 
plastic mounting material with the epoxy, as previously 
described. 

Although there were other varied types of damage, 
the major fracture mechanism was the failure of the zero 
degree lamina oa the side opposite the load. It is 
conjectured that the "field" point encountered is the 
failure of the outer layer in the transverse direction and 
the shifting of this load to tne inner laminas. Due to the 
geometry of the lay-up and the strength characteristics, the 
next two 45-degree laminas could more efficiently carry the 
load transverse to the zero-degree layer. 

The useful, or safe working, load of the plate would 
be limited by this "yield" point. 
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Figure 32 



LIGHT aiGRDSCDPE, "II3LD" 
SEGTION 5B, 75X 



POINT EXAfllNAriDN, 
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Figure 33 



LIGHT KIC333CDPE, "HELD" 
3EGTI0N 5B, 75X 



POIiiT EX AMINATID.'}, 
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Figure 34 



LIGHT JirCRDSCDPS, "YISLD" 
5E:rrON 5B, 23DX 



POINT EXAMINATION, 
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Figure 35 



LIGHT MIC33SCOPE, "YIELD" 
5EGTI0N 5B, 20DX 



POINT EXAHINATIOM, 
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Figure 36 



LIGHT .'lIC.RDSGDPS, "HELD" 
5EZTI0N 5B, 23DX 



POINT EXAMINATION, 



30 



3- M icr oscooig i.x^5kQI.£i2!L 21. Initial 2Z.ll.2i.£ Da aa qe 

Zone 

Early in the impaot loading history, a '’aruise" was 
noted in the imaediata inaaot area. This abnormality coaid 
first be detected using the liguid crystal inspection 
method. As each test progressed, the “bruise" could be 
detected by carefully signting across the face o£ the panel 
where it appeared as a slight bulge. This damaged area 
never grew larger than the immediate area of impact. 

To examine this phenomenon further, a section (1A2| 
was cut from a panel as shown in Fig. 13. Figure 37 is in 
an area just outside the ’bruise" and shows little damage. 
Figure 38 is in tne edge of the pre-failure damage area, and 
has all plies, with the exception of the middle two, 
starting to delaminate. Figure 39 is well into the impact 
area and delamination is extensive. 

The "bruise" dalamination and the typical failure 
differ. In the "bruise" area, more lamina are involved 
extensively without as many associated matrix failures in 
adjacent lamina. There are also no failures of the inner 
two laminas in the transverse direction. 

This damage mechaiism is theorized to be a form of 
spalling caused by the intensive compressive stress wave 
traveling through the thictness of the material on impact; 
Ref. 10 applies. 
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Figure 37 



LIGHT HIC8D3COPE, LOCALIZED DYNAMIC 
SEZriON 1A2, 100X 



DAHAGE, 
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Figure 38 



LIGHT .'HC33S30PE, LDGALIZED DYNAHEZ DAMAGE, 
SECriDH 1A2 
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Figure 39 



LIGHT arcaOSCOPE, localized 
SECTION 1A2, IDOX 



OyNAMIC 



DAMAGE, 



84 



MiS£2sco£ic Exanination of Fiber Failures 



Very few finer failures ware found, and tie cause of 
those failures was hard to determine. 

Figure 40 shows a single bcaicen fiber in the 45 
degree layer of the Scanning Electron Microscape sample. 
Figure 41 shows a bundle af bro>cen fibers in the zero degree 
outer layer of the same specimen. 

The broken fibers in these pictures show 45 degree 
fracture surfaces. This would be indicati/a of shear 
fracture under high strain rates according to Soangelo and 
Heiser [Ref. 11J. Ihef could also just be the ends of the 
fibers as manufactured. 

Because af the rarity af broken fibers and the fact 
that fibers are not normally continuous in a Lamina, no 
significance was assigned to these fiber fractures. 

It appears the lack of strength of the Laminas in 
the transverse direction was the significant factor in 
failure of these plates. Adding 9D degree laminas might 
have provided a large increase in impact resistance and 
static strength. 
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Figara 40 



5EM, A BROKEN FIBER, 520DX 
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BROKEN 



FIBERS, 2200X 



Figure 41 



- SEM, i BUNDLE OF 
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ENERGY REQDIRED FOR FilLURE 



Figure 42 is a plot of the energy rsguired for ultimate 
failure for both dynamic and static cases. Only simply 
supported plates, using the small penetrator, wece plotted. 
A least- squares fit of first and second order were plotted, 
and, within experimental accuracy, the energy for failure in 
both static and dynamic cases was constant. 

Since a constant impact energy mechanism was manifest 
and this energy equalled the static strain energy to 
failure, the correct dynamic model seems to be a one 
degree-of-f reedom lumped parameter spring mass system [Ref. 
9]. The general, continuous, structural dynamic response 
model was not indicated. However, an extensive series of 
dynamic response studies was performed using the dynamic 
response section of SA? E/ [Ref. 5]. The parameters of 
impact duration, impact shape, number of natural modes 
utilized, and constitutive properties were all varied. In 
all results the experimental values showed the actual plate 
as being significantly stiffer than the model, and no 
correlation between the experimental results and the finite 
element dynamic structural response analysis could oe 
obtained . 
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Figure 42 - IMPACT ENE.4GY KEQUIRED FOR FAILURE 
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D. 



STATIC SOLUriDNS 



With a const aat impact anergy failure, a static solution 
to the plate problem would also provide tha dynamic 
solution. Figure 43 is a plot of several different 
solutions and the average experimental results. 

The exact solution of Timoshenko and Woinowsky-Krieger 
■Ref. 7] for an orthotropic plate were computed aid plotted. 
The design manual's [Ref. 2] exact solution was also plotted 
using both sets of developed constitutive properties. 

The finite element model was used for the fixed edge 
solution, and the simply supported solution using the small 
penetrator. It was also used foe the simply supported 
solution using the large penetrator. The large penetrator 
was simulated using the experimentally obtained 
displacements. A drawing was made of the displacements with 
the penetrator overlaid. The load was assumed to be a 
pressure distribution in the same shape and proportion to 
the contact area of the penetrator on the displaced plate. 

As can be seen from Fig. 43, the theoretical results 
agreed well among themselves. The exact and finite element 
solution for a concentrated load were within a percentage 
point; however, the theoretical and experimental results 
were in significant disagreement. 

From the microscopic analysis, the static theoretical 
analysis, and the dynamic theoretical analysis it was 
determined that small deflection theory was unusable. These 
plates, due to their small thickness and relative 
flexability, were more appropriately considered in large 
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deflection theocy 



Chi-Teh Wang [ Ref. 
a large defle:;tLon, 
plate. Although not 
considered here, his re 
a uniformly loaded plat 
concentrated load plate 



12] has solve 
sLnply suppo 
directly ap 
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3* trtU3^- 
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Figure 43 EXACT, FINIfE ELEMENT, AND EXPERIMENTAL RESULTS 

DF TUE STATIC TESTS 



VIEI. 



CONCLUSIDNS 



k. The dynamic ani static failare ncies for simply supported 
graphite-epoxy plates under either central concentrated load 
or low-energy impact are the same. 

B. Between zero and twenty feet per second large deflection 
failures of graphite-epoicy , simply supported plates occurred 
at a constant energy. 

C. There was a "yield point" where damage occurred in static 
loading. 

E. There was localized stress wave damage in the immediate, 
dynamic contact area. 

F. Linear, small-deflection theory is not applicao le- to thin 
graphite-epoxy plates loaded to failure. 
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IX. 



3EC3MMENDMI3N5 



An exact, large deflection solution should be developed 
to analyze thin plate responses to concentrated loads. Once 
this solution is in hand, a good approxination caa be made 
for low energy, dynamic damage thresholds using energy 
methods. This information can then be used for both design 
and operational maintenance. 
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